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ABSTRACT 


A single-stage transonic research compressor and test rig are to be used to obtain 
data on the effect of inlet flow distortion on compressor (and therefore engine) stall. 
Auxiliary injection was examined as a technique for generating distortion in inlet 
stagnation pressure, or temperature, or to simulate the more complex effects of engine 
steam ingestion from a catapult launch. Engineering analyses were developed and 
programmed in EES to relate inlet conditions to the compressor characteristics, for both 
pressure and temperature distortion. An injection duct area of 8% was selected to limit 
the required heater power. A CED analysis was carried out to predict the compressor inlet 
flow field and hence position the injection duct exit. It was found that a broad range of 
distortion parameters could be generated by simply ducting (and heating) atmospheric air 
(or steam) through an auxiliary inlet throttle valve. 
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I. INTRODUCTION 


Engines for military fighter aircraft must be designed to operate stably over a 
required flight envelope. An adequate “stall margin” is usually an engine design 
requirement. Since distortion of the flow into the fan or compressor is known to reduce 
the stall margin, stable operation with a specified level of inlet distortion (due to 
imperfect flow through the inlet) is also usually specified as a design requirement. 

The joint strike fighter (JSF) faces two other potential inlet distortion problems, - 
gas re-ingestion due to jet deflection, and steam ingestion on carrier take-off. The 
potential problem of gas re-ingestion of the STOVE version in hover was lessened by the 
selection of the lift-fan over the jet-lift concept, but the potential for an engine stall on 
carrier take-off must be addressed in view of the single-engine aircraft design. In an 
effort to understand the temperature variation due to steam ingestion during the catapult 
launch, the JSE program office recently conducted a series of tests to survey the 
temperature field and the location of steam exiting a catapult during launch (Donelson S., 
Briggs T., 2003). However, at the current time, steam ingestion from an aircraft carrier’s 
catapult is one of the least understood potential mechanisms for engine stall. 

Most concentration has been placed previously on pressure distortion (produced 
by inlet separations), and temperature distortion (from engine gas re-ingestion). The 
Society of Automotive Engineers (SAE) addressed the aircraft engine-airframe 
compatibility development process, how it was affected by flow distortion, and more 
particularly how to analyze the inlet total pressure distortion effects on gas turbine 
engines (SAE - AIR 1419, 1983). Also, the SAE S-16 committee reported tests to 
examine temperature distortion effects on the intake/engine aerodynamic compatibility, 
and their impact on the design and development of aircraft propulsion systems. Engine 
performance degradation, including power loss due to the compressor instability, has 
been attributed to engine inlet total temperature distortion (SAE ARD50015, 1991). 

The overall goal of the present study is to help develop and validate methods that 
can predict the effects of inlet flow distortion on the fan and compressor, and therefore 
engine, stability. It is clear that the analysis of the aerodynamic response of transonic 
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compressors and engines to non-uniform inlet flow is extremely complex. Consequently, 
if successful methods are to be developed, it is necessary also to conduct controlled 
experiments with distortion, to generate and measure controlled non-uniformity in both 
pressure and temperature into the fan or compressor, and obtain validation data for the 
prediction of instability and stall. 

The need for detailed validation data motivated the present effort, which was to 
design an inlet distortion generation arrangement for the Naval Postgraduate School’s 
transonic compressor test rig. This unique rig, incorporating a throttled inlet flow, 
allowed a somewhat unusual approach to be taken to the generation of distortion, one that 
would allow control of the degree of distortion, whether it was due to air stagnation 
pressure or temperature, or due to steam. The basic idea was to introduce a secondary 
flow with higher stagnation temperature, higher or lower stagnation pressure, or 
containing steam, into the flow entering the compressor. The design process required 
identifying a conceptual arrangement for the generation, developing an analysis of 
parallel-flow compressible mixing, which allowed the specification of the fully mixed- 
out condition, and the computational prediction of the mixing flow after the size of the 
injection duct had been selected. 

In documenting the design, the compressor rig and distortion generation concept 
are first described in Chapter II. The analysis of compressible mixing, given the ratio of 
stagnation pressures, or temperatures, in parallel streams, is then described in Chapter III. 
The unusual requirement here was to be able to specify a required value for the corrected 
flow rate and stagnation pressure after mixing (supplied to the test compressor). The 
analysis allowed the secondary flow duct area to be selected such that it could provide 
pressure or temperature distortion with acceptable auxiliary piping and heater power 
levels. 

The results of analyzing the proposed arrangement using the CFD code ACE are 
then described and discussed in Chapter IV. Conclusions and recommendations are given 
in Chapter V. 
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II. COMPRESSOR AND INLET GEOMETRY 


A. FACILITY AND OPERATION 

A schematic diagram showing the Transonic Compressor Rig (TCR) in the high¬ 
speed building at the Turbopropulsion Laboratory (TPL) at NPS is shown in Figure 1 (O’ 
Brien J.M., 2000). The test compressor is driven by two opposed rotor air turbine stages, 
supplied by a 12-stage Allis-Chalmers (AC) axial compressor. The AC compressor can 
supply air continuously at up to 300 KPa at flow rates up to nearly 5 Kilograms per 
second. Pressurized air from the compressor is fed through a motor driven valve into the 
turbine drive unit. The electric drive motor is manually adjusted to control the 
compressor speed. 

The test compressor draws atmospheric air through a rotating-plate throttle valve 
into a one-meter diameter settling chamber. A 0.46-meter diameter pipe, 5 meters in 
length and containing a flow nozzle, connects the chamber to the test compressor. A 
smooth contraction is provided between the 0.46-meter inlet pipe and the 0.279-meter 
diameter compressor case wall. The flow enters axially into the test compressor rotor and 
then, after the stator, exits axially through a honeycomb flow straightener, as shown in 
Figure 2. Not shown in these figures are an additional 1 MPa compressor that supplies air 
to a balance piston located on the drive shaft, which controls the axial force on the 
bearings of the rotor; also, a smaller shop compressor that provides dry air for the bearing 
oil-mist lubrication system, and other instrument air requirements. 

The Transonic Compressor Rig is operated under manual control from a protected 
control room outside the compressor test cell. Two electrically driven butterfly valves are 
used to adjust the speed of the turbine, while maintaining adequate flow through the AC 
compressor. The rotating plate throttle, which is hydraulically actuated using a solenoid 
valve, is used to set the desired flow rate through the test compressor. The air pressure to 
the balance piston is adjusted by adjusting the instrument air pressure to a pneumatic 
regulator. Hence the operator must adjust four different valves in order to set the 
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operating condition for a single point on a constant 
This is significant when the proposal is to add aux 
additional valve. 



Figure 1. Schematic diagram of the transom 


TEST COMPRESSOR 










B. TRANSONIC COMPRESSOR 

The single stage transonic compressor was designed by Nelson L. Sanger at the 
NASA Glenn Research Center using CFD techniques (Sanger, 1996). The highly loaded 
design resulted in 22 blades in the rotor and 27 in the stator. The performance of the stage 
was mapped experimentally in the TCR (Gannon et al, 2004), and predicted using CFD 
codes (Hobson et al, 2004). Fixed instrumentation for performance mapping was installed 
at the three stations shown in Figure 2, and in the 0.279-meter diameter duct ahead of the 
rotor. It is significant to the present work that the case wall can be rotated to survey 
peripherally using fixed probes. 

C. COMPRESSOR PERFORMANCE MAP 

A performance test is usually conducted while maintaining a constant compressor 
speed. The test begins with the main throttle fully open. Data are then taken at various 
main throttle positions, giving desired increments in the flow rate. The stall boundary is 
established at low rotational speeds and is approached cautiously at near-design speeds. 
When stall does occur, it is necessary to open the throttle to restore stable operation, and 
possibly close it gradually again to more closely approach the stall boundary. 

Performance map data are given Appendix A. (from Gannon, 2004) and data for 
80 %, 90% and 100% of design speed are shown plotted in Figure 3. The data in Figure 3 
represent the pumping characteristic of the installed Sanger stage when operating without 
distortion. It is likely that less pressure ratio will be produced at a given corrected flow 
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when distortion is introduced. This was considered in designing the auxiliary injection 
system to produce controlled distortion. The 90% speed line was used rather than the 
100% speed line, and a 5% degradation in the pumped pressure was assumed at the 
compressor inlet. 




* 80% Speed 
eperimental data 


• 90% speed 
experimental 
data 


♦ 100% speed 
experimental 
data 


Figure 3, Experimental compressor performance (without inlet distortion) 


D. DISTORTION GENERATION CONCEPT 

Auxiliary injection into the inlet pipe, in order to generate controlled magnitudes 
of well-defined inlet distortion, is proposed. In order to avoid three-dimensional effects 
that would both complicate the prediction and require extensive inlet flow diagnostics, 
injection parallel with the main flow is proposed. A sketch of the conceptual arrangement 
is shown in Figure 4. 
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Section A-A 




b) Duct detail 

Figure 4, Auxiliary injection concept for generating controlled inlet distortion 
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III. ANALYSIS OF AUXILIARY INJECTION 


A. GEOMETRY AND NOTATION 

Figure 5 shows the geometry that was used. Station 1 represents the main flow; 
station 2 represents the secondary or distorting flow; station 3 represents the end of the 
mixing of the main and secondary flows, and station 4 represents the entrance of the 
compressor. 



Figure 5. Schematic representation of the distortion generation hy secondary flow 

injection 


Notation: 

M Mach number 

m Mass flow rate 

P Pressure 

Pt Stagnation pressure 

T Temperature 

Tt Stagnation temperature 

V Velocity 
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R Gas constant 

Cp Specific heat at constant pressure 

Cj, Specific heat at constant volume 

y Ratio of specific heats (here, y = lA) 

a Ratio of secondary to overall flow area (a = —) 

Aj 

Pt 

71 Ratio of secondary to main-flow stagnation pressure {n = —^) 

Pt^ 

Tt 

T Ratio of secondary to main-flow stagnation temperature {t = —-) 

Tt^ 

Pt^ 

V Ratio of stagnation pressure at station 3 to design value (V = — ^) 

Suffixes: 

1 Main flow 

2 Injected flow 

3 Fully mixed-out flow (in inlet pipe) 

4 Compressor inlet flow (after bell-mouth contraction) 

d Design condition (no distortion) 

std Standard atmosphere 

B. PRESSURE DISTORTION GENERATION 

1. Analysis 

The following assumptions were made: 

1. Neglect effect of friction at the wall 

2. The two flows are at the same stagnation temperature (x = 1) 

3. After mixing, the flow is uniform at station 3. 

4. The flow from station 3 to station 4 is isentropic ( Pt^ = Pt^ , Tt^ = Tt^) 

5. The injected stream enters parallel to the main stream, so that Pi = P 2 . 
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It is assumed that Pt 3 is a fraction of Pt 3 d, and that the fraction can be greater or 
less than unity (positive or negative pressure distortion). The goal is then to find the 
magnitude of Pt 2 which results in the required value of Pt^. 

Using conservation of mass 

m3 = m2 + mj (1) 

From conservation of energy, in the absence of shaft work and heat addition, 
assuming a perfect gas, 

m^Tt^ =m2Tt2 +m^Tt^ (2) 

From conservation of momentum, 

0 = mjVj + mjV 2“^3^3 + ^3 [^1 “ A ] (3) 

The area ratio a , is given by, a = — (4) 

4 

and a value 0< cr <1 is selected. The magnitude of the pressure distortion is defined as 
Pt 

—- = TT . Both .;r>l and .;r<l are of interest. 

Pt, 

The analytical approach followed was to express all the quantities with respect to 
the unknown Mach number at station 1 {M,). The development is described as a series of 
steps: 

1. Assume a value of Pt^lPt^^ • The value of Pt^,^ is taken from experimental data 
presented in Appendix A for 90% of the design of speed. The ratio is defined as 

= V , (V<1 or >1) (6) 

2. Evaluate Tt^ from conservation of energy (equation 2) to obtain 
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(7) 


Ttj^ 


m^Tti + 771^^2 


m. 


Since= Tt^, using equation (1), 

rtj = Tt^ = Ttj, 


( 8 ) 


3. At any station, the flow rate can be expressed in terms of stagnation conditions 
and Mach numbers as in Appendix B. 


m = 


PtA 

{yRTtf^ 




1 + 




(r+i) 

2(r-l) 


(9) 


Pt^ 


Using the conventional (engine) notation, d = —7 ^ where =standard 


std 


Pt 

atmosphere, ^ is obtained using equation (9) and equation (1) as 

Pt^ 


J31 — 7 t\ 


M, 








ir+i) 

^ 2(7-1) A. 

+ —M, 


) 




M, 


(7+1) 

2(7-1) 


_ (7+1) 

2(7-1) 


M. 


l^(7zl) 




( 10 ) 


4. While 1^ higher or lower thanPtj , since both flows are subsonic P 2 
Therefore Mj is different from M 2 , but the two Mach numbers are related by 


1 + 



7 

1 + 

V 



) ¥ 
\n ^ 


( 11 ) 


5. Finally, using equations (9) and (10) in equation (3) and introducing the 
relationships between velocity and Mach number, static and stagnation temperature, the 
following is obtained from the conservation of momentum. 
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o = s,,r\ 




(r+i) 


\ + K - ^ 


A, M/ 


(y+O n, 

2 + 1)m/^ Ax-i) 


2+(r-l)M,^ 








(y+i) 

'2(r-l) 




M 


7 

2 


( 12 ) 


where (^,, = — 

' Pt^ 

Equation (11) and equation (12) can be written as a single equation for one 
unknown Mach number. An EES code was written to solve first for Mj and then M 2 and 
then calculate all the other unknowns (Appendix C). 

2. Results 

The inlet distortion design calculations were based on the assumption that 

Pt. 

— = V = 1.05 (that the compressor pumping capability would deteriorate with 


flL . ± I ± , 

distortion present), and the corrected flow rate, ^-—, would be given by the 

! ^std 

experimentally determined compressor map shown in Eigure 3. 


At 90% of the compressor design speed, the ratio of the inlet pressure to standard 
reference pressure measured in tests, with a second curve representing 5% deterioration, 
are shown in Eigure 6. The unusual trend in the behavior of inlet pressure with throttling 
is the result of the throttle being in the inlet flow to the compressor. The upper curve was 
taken as the design curve for the distortion generation. Similar curves were calculated for 
80% and 100% speed levels. 
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Figure 6, Compressor pumping characteristic measured at 90% design speed and 

degradation assumed for distortion design. 

For reasons given in the following section, a design value of q;= 0.08 was 
selected. Calculations were then carried out for pressure distortion values of n from 
slightly smaller than unity to 1.05. Since the Mach number in the inlet pipe is very low, 
the ratio of the static to stagnation pressure (Pi/Pti) is very close to unity, and the 
minimum value of 7 i (zero injection velocity) is, from Equation 7 of Appendix E, 
71=0.9965. Hence the full range of ‘negative pressure distortion’ is achieved by varying n 
only a fraction of a percent below unity. Similarly, the complete range of interest for 
‘positive pressure distortion’ is achieved in the range l<7l<1.05. Since little information 
is lost by the omission, only the results for 7i > 1 are shown plotted in the results. 

The results are given in Eigures 7-9 for operation at 80%, 90% and 100% of 
design speed. The Mach numbers in the main and secondary streams are shown in Eigure 
7, the stagnation pressures at stations 1 and 2, required to satisfy P 73 =1.05 are 
shown in Eigure 8 , and the mass flow rates in the two streams are shown in Eigure 9. 
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Note that all the blue curves are for the near-stall corrected mass flow at 90% operational 



Figure 7, Mach number variation with pressure distortion parameter 

(a =0.08, v = 1.05) 



Figure 8, Injection stagnation pressure variation with distortion parameter it 

(« =0.08,1/= 1.05) 
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- 80% speed 


90% speed 


100 % speed 


-Total mass flow 
rate 80% speed 


-Total mass flow 
rate 90% speed 


-Total mass flow 
rate 100% 
speed 


Figure 9, Mass flow-rate variation with distortion parameter (a =0,08, v=l,05) 


Figure 7 shows that when ;r increases, M^ rises at a much higher rate than the rate 
at whichMj decreases. Similarly, in Figure 9 increases and decreases with 
increasing distortion. The required stagnation pressure in the injected flow Ptj > increases 
with 71 but remains lower than standard atmospheric pressure (101325 Pa). Therefore, 
the injected flow can be fed from the atmosphere and the pressure controlled using a 
throttle. With the throttle in the injected flow, the amplitude of Pt^ can be adjusted, 
while the throttle in the main flow can be used to adjust the overall level of Pt ^. 
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c. 


TEMPERATURE DISTORTION GENERATION 


1. Analysis 

The following assumptions were made: 

1. Neglect effect of friction at the wall 

2. The two flows are at the same stagnation pressure (7l = 1) 

3. After mixing, the flow is uniform at station 3. 

4. The flow from station 3 to station 4 is isentropic ( Pt^ = Pt^, Tt^ = Tt^) 

5. The injected stream enters parallel to the main stream, so that Pi = P 2 . 

With these assumptions, the equations for conservation of mass, momentum and 

energy remain the same as given for the analysis of pressure distortion in Equations 1, 2 
and 3 respectively. Also, as before, the area fraction a in Equation 4 must be selected to 
have a value 0<C1;<1. Eurthermore, it is again assumed that Pt 3 is a fraction of Pt 3 d 
Pt, 

(—^ = v >1) to account for an expected degradation in pumping capability with inlet 

distortion. Einally, since the main flow will always be pumped from the atmosphere, it 
will be assumed thatTtj = . 

Tt 

The magnitude of the temperature distortion is defined by the value of T = —-. 

Tt^ 

While both x >1 and t< 1 are of interest, x >1 can be implemented more easily, using a 
heater in the auxiliary flow. Since the stagnation pressures are the same, and the static 
pressure is common, the Mach number in the main flow is the same as in the injected 
flow. Thus Ml = M 2 , T 2 = T^t , and the entering velocities are different. 

The analytical approach followed was again to eventually express all the 
quantities in terms of the unknown Mach number at station 1, M^. The development is 
described in a series of steps. 

1. Assume a Ptj smaller or bigger than Pt^^. The value of is taken from 

experimental data presented in Appendix A for 90% operational level of speed. 
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= V <1 or >1 


(13) 




Note that the corrected flow rate is then known from the compressor map. 


2. Using Equation 9 at stations 1 and 2, with 7 l = 1, Mi = M 2 and T given, and writing 

(1-q;) 

a 

mj ={X4t)m^ (14) 


and using Equation 1, 




^ Ayfr + 1 ^ 




(15) 


3. Using Equations 14 and 15 in Equation 2, 

Vt, 


Tt^ = 




1 + A \ T 


(16) 


Since Tta is now known, using the corrected flow rate and Pta, the flow rate m 3 is 
obtained; then mi and m 2 are given by Equations 15 and 16, respectively. 


4, Using Equation 9 (at stations 1 and 2) in Equation 3, 


f ^ 

— + A3 

{RTtJ^ 

V 

1+ ^ M, 

1 1 2 J ‘ J 

0.5 

/ 

M,\ 

( 2) 

1+ ^ M, 

1 1 2 j ‘ J 

{h+b) 3 


(17) 


where the exponents are h = 


-ir+i) 

2(r-i) 


and 


b = 


7 

r-i' 


Since V 3 and P 3 in Equation 17 can be written in terms of stagnation pressure and 
temperature, which are known, and Mach number, which is obtained from Equation 9, 
the only unknown is Mi. 
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An EES code was written to evaluate first Mj and then calculate all the other 
unknowns (Appendix D). 

2. Results 

The inlet temperature distortion design calculations were again based on the 
Pt. 

assumption that -= V = 1.05 (that the compressor pumping capability would 

deteriorate with distortion present), and that the corrected flow rate, m^ = 

would be given by the experimentally determined compressor map shown in Eigure 3. 
Erom preliminary calculations of the power required for different areas of the injected 
flow, a value of «=0.08 was selected. Results were then obtained for temperature 
distortion values of Tfrom 1.05 to 1.25.The results are given in Eigures 10-12 for 
operation (near stall) at 80%, 90% and 100% of design speed. 

The Mach number in the main and secondary streams is shown in Eigure 10. The 
amount of energy transferred to the injected flow at stations 1 and 2, required to satisfy 
Pt^ = \.Q5 Pt^^, is shown in Eigure 11, and the mass flow rates of the two streams are 

shown in Eigure 12. Note that all the blue colored curves are for the minimum design 
corrected mass flow rate at 90% operational speed. 


PJP.d 
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Figure 10, Mach number variation with temperature distortion parameter t {a =0,08, 

v=l,05) 
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Figure 11, Energy required to heat the injected flow (« =0,08, v=l,05) 
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Figure 12, Mass flow-rate variation with distortion parameter x{a =0,08, v=l,05) 


Figure 10 shows that when the magnitude of r increases, Mj (and therefore M 2 ) 
remains approximately constant. Similarly, as x increases, there is little effect on the mass 
flow rates in the two streams. This is to be expected since it has been assumed that the 
stagnation temperature in the main flow will not change. The energy required to heat the 
injected flow is seen to increase directly withr. The data in Figure 11 explains the 
particular selection made for the injection duct area ratio. With the selection of a=0.08, it 
was possible to vary x up to 1.25 without exceeding 25 kW in the required power. 
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IV. COMPUTATIONAL ANALYSIS OF DISTORTION 

GENERATION 


In the previous chapter the analysis of the auxiliary injection assumed that 
uniform (“fully mixed out flow”) conditions occurred at the downstream station, and this 
allowed the downstream boundary condition to be related to the compressor-pumping 
characteristic (or “map”) at Station 4. The analysis properly conserved mass, momentum 
and energy, and served to allow the size of the auxiliary duct exit to be selected as shown 
in Figure 13. In reality, the auxiliary and main flows will not be fully mixed out, and the 
flow at the compressor face will depend on the position of the auxiliary duct exit forward 
of the compressor; i.e. on the selection of the “mixing length”, L, shown in Figure 14. 


28 . 8 ° 
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A computational analysis of the compressible mixing was carried out using the 
CFD code ACE (References 8, 9, and 10). Figure 14 shows the geometry that was 
modeled in order to evaluate the mixing for 3 different locations of the exit of the 
auxiliary duct; specifically, for L=0.5, 1.0 and 2.0 meters. 


Contraction 



Figure 14. Geometry of the duct used in CFD ACE code, 

ACE is one of a set of computer codes for multi-physics computational analysis 
developed by CFDRC. The codes provide an integrated geometry and grid generation 
module, CFD-GEOM; a graphical user interface for preparation of the model, a 
computational solver for performing the simulation, CFD-ACE (U), and an interactive 
visualization program for examination and analysis of the simulation results, CFD- 
VIEW. 
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There are provisions to create a computational grid with boundaries that can be 
classified as symmetry, inlet, outlet, wall or ‘arbitrary’. Stagnation pressure and 
temperature can be specified and held at the inlet, static pressure can be specified and 
held at the outlet. In order to compute the flow corresponding to 7l = 1.05 (pressure 
distortion) or T = 1.25 (temperature distortion), the inlet stagnation pressure level was 
adjusted until the mass flow rate was equal to the required value. 

An example of the computational grid generated to model the main duct with an 
auxiliary inlet duct occupying 8% of the annulus is shown in Figure 15. Note that a 
rectangular region was required along the axis to avoid the singularity from a center axis. 



Figure 15, Example of the grid used to calculate the mixing from auxiliary injection 

When the grid was satisfactory, a converged solution (three orders of magnitude 
reduction of residuals) could be obtained, as shown in Figure 16 for one case of pressure 
distortion. Figure 16 shows the larger inlet pipe, with uniform inflow over the auxiliary 
injection sector, and the mixed profile at the exit (compressor face) after the area 
contraction. In order to compare the profiles generated with different mixing lengths, data 
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at the outlet boundary were plotted along the two lines shown in Figure 16 and 
dimensioned in Figure 17. Results obtained for pressure and temperature distortion are 
described in the following sections. 


Measurements on 
Z axis z=0.07m 


P_tot - N/m^2 

8.391 E+004 



8.3E+004-« 

8.2E+004- 

8 1E+004- 

I 

BE+004-^ 

7.9E+004- 

7.8E+004- 

7.776E+004 


Figure 16, Inlet to exit flow and lines chosen at the exit for data comparison 



Figure 17, Locations that were selected at the exit plane (station 4) to compare results 
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A. PRESSURE DISTORTION 

The results of calculating the mixing with the auxiliary flow injected at 5% higher 
in stagnation pressure (n = 1.05), but at the same stagnation temperature (x =1), as the 
main flow, for three different mixing lengths, are shown in Figures 18 and 19. 


Figure 16 and 17 shows the locations that were selected at the exit plane (station 
4) to compare results. It is observed that the levels of inlet stagnation pressures were 
higher than these that were used in corresponding cases in chapter III, since the exit static 
pressure and flow rate were the same but the exit was not fully mixed out. 
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Figure 18, Stagnation pressure distribution at station 4 (at Z = 0,044 m) for different 

inlet duct lengths (90% speed, n = 1,05, x = 1) 
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Figure 19, Stagnation pressure distribution at station 4 (at Y = 0,115 m) for different 

inlet duct lengths (90% speed, n = 1,05, x = 1) 


B, TEMPERATURE DISTORTION 

The results of calculating the mixing with the auxiliary flow injected at 25% 
higher stagnation temperature (x = 1.25), but at the same stagnation pressure (n = 1.0) as 
the main flow, for three different mixing lengths, are shown in Figures 19 and 20. Figure 
16 shows the locations that were selected at the exit plane (station 4) to compare results. 
Again it is observed that the levels of inlet stagnation pressures were higher than those 
that were used in corresponding cases in chapter III, since the exit static pressure and 
flow rate were the same but the exit was not fully mixed out. 
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Figure 20, Stagnation temperature distribution at station 4 (at Z = 0,044 m) for 
different inlet duct lengths (90% speed, 7t = 1,0, x = 1,25) 
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Figure 21, Stagnation temperature distribution at station 4 (at Y = 0,115 m) for 
different inlet duct lengths (90% speed, 7t = 1,0, x = 1,25) 
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From Figures 18-21, it can be seen that the mixing in the case of temperature 
distortion is less than for pressure distortion, although the magnitude of the distortion 
parameter is larger (x = 1.25 compared to 7l = 1.05). In fact, what drives the mixing is the 
difference in the incoming velocities of the two streams. It can be shown (Appendix E) 


that 



when n = 1.0, but that — 



when T = 1.0. 


Since the Mach number in the inlet pipe is low (~ 0.07), 


= 3.85 for the pressure 


distortion case analyzed (n = 1.05) compared to 1.12 for the temperature distortion case 
analyzed (x = 1.25). These values of relative inlet velocity were confirmed by the code. 


C. GRID SENSITIVITY 

The sensitivity of the CFD code predictions to different grid selections was 
evaluated for the 0.5 meter mixing length geometry. Three different grids were used, 
progressively increasing the number of nodes. Grid 1 had 37,908 nodes. Grid 2 had 
85,183 nodes and Grid 3 had 128,478 nodes. The results are shown in Figures 22-25. It 
can be seen (in Figure 23) that the largest effect of the grid selection was from Grid 1 to 
Grid 2 in the case of pressure distortion. However, from Grid 2 to Grid 3 there were very 
small differences in the results for both pressure and temperature distortion. 
Consequently, Grid 2 was used for the 0.5 meter mixing length, and a similar grid density 
was used when extending the inlet duct length (100,103 nodes for the one meter and 
133,366 nodes for the two meter length). 
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Figure 22, Stagnation pressure distribution at station 4 (at Z = 0.044 m) for different 

grid selection (90% speed, 7t = 1,05, x = 1) 
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Figure 23, Stagnation pressure distribution at station 4 (at Y = 0,115 m) for different 

grid selection (90% speed, 7t = 1,05, x = 1) 
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Figure 24, Stagnation temperature distribution at station 4 (at Y = 0,115 m) for 
different grid selection (90% speed, 7t = 1,0, x = 1,25) 






Figure 25, Stagnation temperature distribution at station 4 (at Z = 0,044 m) for 
different grid selection (90% speed, it = 1,0, x = 1,25) 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The present study was motivated by the need to conduct controlled experiments on 
the effect of steam or hot-air ingestion on the stability and stall margin of gas-turbine 
engine compressors. The proposal to use auxiliary injection into the inlet pipe of the NFS 
transonic compressor was examined. The following were concluded: 


• It is feasible to generate and control pressure and temperature distortion using air 
ingestion through an auxiliary duct incorporating one additional auxiliary throttle 
valve. 

• Using a 25 KW heater, a suitable range of temperature distortion can be 
developed using an auxiliary duct area equal to 8% of the inlet pipe. 

• Two engineering codes (EES) were generated to calculate inlet conditions from 
the compressor characteristic for specified pressure and temperature distortion 
parameters. The code calculations enabled the auxiliary duct size to be selected. 

• A CED analysis was carried out successfully to determine the flow profile 
entering the compressor and the effect of the available mixing length. The results 
were shown to be insensitive to grid selection. 


The following recommendations are made: 

• Extend the CED simulation to include the full geometry of the injection duct 

• Complete the inlet hardware design (with provision to position the auxiliary duct 
at different locations) 

• Design, or select, a suitable air heater 

• Procure a steam generator to inject through the same duct. 

• Make provision for the measurement of secondary air/steam mass flow rate, 
pressure and temperature. 
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APPENDIX A TEST AND INLET DISTORTION DATA 


A. EXPERIMENTAL PERFORMANCE MAP 


Table 1, 80% operational speed data 


Run no 

1 

Mass Flow 

Pt ratio 

Tt4d(0c) 

Pt4d 

P4d 

RPM 

M4 

2 

3.99164818 

1.344786 

16.8304362 

82257.85 

79766.09 

21827.19 

0.210088 

3 

4.19552573 

1.337583 

17.69664714 

83233.64 

80499.93 

21799.52 

0.218943 

4 

4.583530331 

1.322167 

17.61668294 

85164.86 

82133.92 

21791.87 

0.228117 

5 

4.945168818 

1.318621 

17.71622721 

86862.63 

83228.69 

21786.58 

0.247841 

6 

5.203798788 

1.309887 

17.02320964 

88343.52 

84513.06 

21755.76 

0.252441 

7 

5.814255239 

1.27468 

17.80014648 

93667.14 

89144.06 

21794.26 

0.266848 

8 

6.055224276 

1.258274 

17.17148438 

96307.87 

91439.04 

21793.18 

0.273245 

9 

6.178094964 

1.252228 

17.30166016 

97353.23 

92371.1 

21769.75 

0.274975 


Table 2, 90% operational speed data 


Run no 








1 

Mass Flow 

Pt Ratio 

Tt4d(0c) 

Pt4d 

P4d 

RPM 

M4 

2 

4.381216587 

1.46121522 

29.4043457 

76802.92 

73402.47 

25086.49 

0.255173 

3 

4.468290403 

1.45931496 

14.98787435 

76679.29 

73529.92 

24394.9 

0.245507 

4 

4.512427704 

1.45328577 

15.70483398 

77083.1 

73692.37 

24380.09 

0.25432 

5 

4.839665367 

1.4535476 

29.01258138 

77784.41 

74220.32 

24922.73 

0.259721 

6 

5.080933861 

1.43758367 

16.29047309 

79539.32 

75537.24 

24434.66 

0.272576 

7 

5.345904496 

1.4294289 

28.50551758 

81149.71 

76826.74 

24974.81 

0.280747 

8 

5.661951478 

1.41257914 

14.23828125 

83281.68 

78574.72 

24362.89 

0.289494 

9 

5.681440546 

1.40218943 

28.14306641 

85147.89 

79734.49 

25014.69 

0.307775 

10 

6.034647541 

1.37024027 

29.19163411 

88890.65 

83267.79 

24946.83 

0.306964 

11 

6.036431418 

1.38713393 

27.87753906 

87052.26 

81500.33 

24916.78 

0.308281 

12 

6.220627881 

1.37083823 

15.25431315 

88720.71 

83176.34 

24373.44 

0.305027 

13 

6.697960407 

1.32553477 

29.42128906 

95424.21 

89074.93 

25078.21 

0.315178 

14 

6.873528918 

1.32203335 

14.99145508 

95457.43 

89226.49 

24323.46 

0.312039 
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Table 3, 100% operational speed data 


Run no 








1 

Mass Flow 

Pt Ratio 

Tt4d(0c) 

Pt4d 

P4d 

RPM 

M4 

2 

4.783778964 

1.593625 

17.10299479 

71356.82 

67204.76 

27146.95 

0.293903 

3 

4.878505613 

1.574358 

20.42115885 

72519.86 

68110.02 

27072.63 

0.300718 

4 

4.964007749 

1.57407 

17.38203125 

73093.4 

68717.47 

26899.6 

0.298286 

5 

5.276035803 

1.571927 

15.20699685 

74177.01 

69238.64 

26940.92 

0.315272 

6 

5.293451804 

1.55561 

18.13105469 

74840.54 

70316.04 

26803.71 

0.299806 

7 

5.378261808 

1.554052 

21.20722656 

74850.87 

68825.11 

27103.79 

0.34835 

8 

5.417627178 

1.561808 

16.52239583 

74722.02 

69676.13 

27065.13 

0.317625 

9 

5.629933906 

1.547752 

18.37011719 

76626.1 

71312.47 

27067.99 

0.322068 

10 

6.060280342 

1.509242 

19.17338867 

80292.57 

74308.21 

26936.42 

0.334491 

11 

6.076155446 

1.517941 

15.41518555 

79847.25 

73591.53 

26972 

0.343395 

12 

6.346813703 

1.479708 

20.77011719 

82624.98 

76112.9 

26978.9 

0.344476 

13 

6.370598575 

1.476179 

17.33754883 

83590.99 

76989.97 

27028.21 

0.344829 

14 

6.48758399 

1.481319 

18.94223633 

83939.83 

77282.39 

27102.85 

0.34562 

15 

6.487651075 

1.48281 

15.18090123 

83668.5 

76948.11 

26965.11 

0.347934 

16 

6.749165973 

1.455565 

15.13730469 

86549.18 

79670.71 

26983.19 

0.345994 

17 

7.06954336 

1.415859 

18.0215332 

90534.76 

83122.17 

27132.38 

0.351473 

18 

7.262751202 

1.388472 

21.27998047 

93822.34 

86224.78 

27121.16 

0.349433 

19 

7.278570267 

1.399183 

18.88925781 

93361.07 

85801.28 

27002.95 

0.349423 

20 

7.400624968 

1.389324 

18.53776042 

94515.02 

86937.39 

26934.37 

0.347596 

21 

7.458758542 

1.393771 

15.44352214 

94426.45 

86798.42 

27020.12 

0.348987 


B. PRESSURE DISTORTION 

In the following tables small changes to the used so far symbols are made, and are 
presented from the notation below. 


Notation: 


a 

P 

mt 


Ratio of secondary to overall flow area (a= 


A2 


) 


Ratio of secondary to main-flow stagnation pressure {p = 
Total mass flow rate mt=m 1 -rm2 


Pt2 

PtT 


) 
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Table 4, Mass flow vs, pressure ratio p 


p 

a 

ml (80) 

ml (90) 

m1(100) 

1X2(80) 

rr2(90) 

1X2(100) 

mt(80) 

rxt(90) 

rxt(IOO) 

1 

0.08 

3.856 

4.232 

4.621 

0.3353 

0.368 

0.4018 

4.1913 

46 

5.0228 

1.01 

0.08 

3.556 

3.986 

4.409 

0.635 

0.6138 

0.6143 

4.191 

4.5998 

5.0233 

1.02 

0.08 

3.355 

3.812 

4.251 

0.8358 

0.7885 

07722 

41908 

4.6005 

5.0232 

1.03 

0.08 

3.194 

3.669 

4.12 

0.9973 

0.9313 

0.9033 

41913 

4.6003 

5.0233 

1.04 

0.08 

3.055 

3.545 

4.005 

1.136 

1.055 

1.018 

4.191 

46 

5.023 

1.05 

0.08 

2933 

3.435 

3.903 

1.259 

1.165 

1.12 

4.192 

46 

5.023 


Table 5, Macb number vs pressure ratio p 


P 

a 

Ml (80) 

Ml (90) 

M1(100) 

M2(80) 

M2(90) 

M2(100) 

1 

0.08 

0.07226 

0.08688 

0.1002 

0.07226 

0.08688 

0.1002 

1.01 

0.08 

0.06668 

0.08188 

0.09559 

0.1367 

0.1448 

0.153 

1.02 

0.08 

0.06297 

0.07836 

0.09224 

0.1799 

0.1859 

0.1922 

1.03 

0.08 

0.06001 

0.0755 

0.08948 

0.2146 

0.2194 

0.2247 

1.04 

0.08 

0.05747 

0.07303 

0.08708 

0.2443 

0.2485 

0.253 

1.05 

0.08 

0.05523 

0.07084 

0.08494 

0.2707 

0.2744 

0.2784 


Table 6, Injected flow stagnation pressure vs, pressure ratio p 


P 

a 

P80 

P90 

P100 

1 

0.08 

86371 

80643 

74925 

1.01 

0.08 

87143 

81368 

75601 

1.02 

0.08 

87897 

82077 

76263 

1.03 

0.08 

88642 

82778 

76918 

1.04 

0.08 

89383 

83473 

77567 

1.05 

0.08 

90120 

84165 

78213 


C, TEMPERATURE DISTORTION 


Table 7, Mass flow vs, temperature ratio t 


t 

a 

ml (80) 

ml (90) 

m1(100) 

m2(80) 

m2(90) 

m2(100) 

mt(80) 

mt90) 

mt(IOO) 

1.05 

0.08 

3.856 

4232 

4621 

0.3272 

0.3591 

0.3921 

5.734 

6.396 

6.839 

1.1 

0.08 

3.856 

4232 

4621 

0.3197 

0.3509 

0.3831 

5.723 

6.384 

6.826 

1.15 

0.08 

3.855 

4231 

462 

0.3126 

0.3431 

0.3746 

5.713 

6.373 

6.814 

1.2 

0.08 

3.855 

4231 

462 

0.306 

0.3359 

0.3667 

5.704 

6.362 

6.802 

1.25 

0.08 

3.854 

423 

4619 

0.2998 

0.329 

0.3592 

5.694 

6.351 

6.791 
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Table 8, Mach number vs, temperature ratio t 


t 

a 

Ml (80) 

Ml (90) 

M1(100) 

M2(80) 

M2(90) 

M2(100) 

1.05 

0.08 

0.07225 

0.08688 

0.1002 

0.07225 

0.08688 

0.1002 

1.1 

0.08 

0.07225 

0.08688 

0.1001 

0.07225 

0.08688 

0.1001 

1.15 

0.08 

0.07224 

0.08687 

0.1001 

0.07224 

0.08687 

0.1001 

1.2 

0.08 

0.07223 

0.08686 

0.1001 

0.07223 

0.08686 

0.1001 

1.25 

0.08 

0.07222 

0.08684 

0.1001 

0.07222 

0.08684 

0.1001 


Table 9, Heat transfer to the injected flow vs, temperature ratio t 


t 

a 

Q80 

Q90 

Q100 

1.05 

0.08 

4766 

5457 

5717 

1.1 

0.08 

9311 

10663 

11170 

1.15 

0.08 

13659 

15642 

16385 

1.2 

0.08 

17826 

20414 

21384 

1.25 

0.08 

21829 

24998 

26186 
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APPENDIX B MASS FLOW RATE EQUATION 


The (uniform) mass flow rate at any section with area A is given by 

m = pAV (1) 

Substituting for the density using the equation of state, 

P 


P = 


RT 


( 2 ) 


and for the velocity using the definition of Mach number, 

V=M^^ (3) 

and then introducing the relationships between static and stagnation properties. 


P = Pt 


\ 


I 2 y 


(r-i) 


( 4 ) 


and 


T =Tt 


1 + 

V 2 y 


-1 


( 5 ) 


Equation 1 can be written as 


(r-l) 

m = , l + 2r-f 




t 


V 


2 y 


(r+i) 
^ 2(r-i) 


( 6 ) 
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APPENDIX C PRESSURE DISTORTION (EES CODE) 


In the EES codes listed in Appendices C and D some changes were made to the 
notation used previously. The following lists departures from the previous notation: 

Notation: 


u Velocity 

k Ratio of specific heats (here, k=1.4) 

A2 

a Ratio of secondary to overall flow area (a= —) 

A3 

Ptl 

p Ratio of secondary to main-flow stagnation pressure {p = ) 

Ttl 

t Ratio of secondary to main-flow stagnation temperature {t = ) 

Pt3 

n, p2 Ratio of stagnation pressure at station 3 to design value (n = p2 = -) 

Pt3d 
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m4d 

= 4.381 

Pt4d 

= 76803 

D4 = 

11 

12 ■ 

3.28 

A4 = 

71 ■ ■ 

D4^ 

4 

K = 

1.4 


R = 

287 


cp = 

R ■ 

K 

K - 1 

TT = 

29.4 


Tt4d 

= TT 

+ 273.15 

D3 = 

18 

12 ■ 

3.28 

A3 = 

D3^ 


4 

“ ‘ 71 

a = 

0.08 



A2 = 

A2 = 

Tt1 = 
Tt2 = 

Tt3 = 


a ■ AS 

A1 ■ 

Tt4d 

Tt1 


1 


Q • 


d 

ll - aj* ' 


1 


1 - a 

1 + - 

a ■ p 


■ Tt1 
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Pt2 = p ■ Pti 
p = Press ureratio 


h = 

-(K + 1 ) 

2 ■ (K - 1) 

m4cl 

M 

= Pt4d ■ - 


Tt4d 

Pt3d 

= Pt4d 

Tt3d 

= Tt4d 

m3d 

= m4d 

m3d 

A3 

= Pt3d ■ - 


Tt3d 

Ma3 

= Ma3d 

P2 = 

= variable 

Pt3 

= P2 ■ Pt3d 

m3 

Pt3 

= m3d ■ - 


Pt3d 


0.5 


0.5 


'K ■ 

0.5 

■ 

fK - 1 1 

0 

R 

■ Ma4d ■ 

1 + 

1 2 i 

■ Ma4d 


'K ■ 

0.5 

■ 

fK - 1 1 

0 

R 

■ Ma3d ■ 

1 + 

1 2 j 

■ Ma3d 


K 


K - 1 


c2 
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K 


1 


C3 = pL ^ J - 1 

c4 = K ■ R ■ Tt1 

Ai r K 1°-^ 

Tt1 _ 


A2 r K 



m3 

p ■ A2 ■ 

1 + c1 ■ 

c3 2 

— + c2 ■ Ma1 
c1 

/ 

c3 2 

— + c2 ■ Ma1 

C1 

0.5 ^ 

+ A1 ■ (1 + c1 ■ Ma1 'Mai 
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P1 = 


Pt1 


(1 + c1 ■ Ma1 ) 


h 


m1 = Pt1 ■ c5 ■ Ma1 ■ (1 + c1 ■ Ma1 ) 
m2 = p ■ Pt1 ■ c6 ■ Ma2 ■ (1 + c1 ■ Ma2 

v1 = Ma1 


v2 = Ma2 


P3 = 


v3 = Ma3 


• 

c4 


0.5 

1 + 

cl ■ 

Mai ^ 


r 

c4 

1 

0.5 

1 + 

cl ■ 

Ma2 ^ 


Pt3 




c1 • 

Ma3 ^ 

C 

) 


r 

c4 


0.5 

1 + 

cl ■ 

Ma3 ^ 


c5 ■ 

Mai 

■ (1 

2 ^ 

+ cl ■ Mai ) 

41 ■ 

c6 ■ 

Ma2 ■ 

( 1 + cl ■ Ma2 


Ma1 


c4 


nO.5 


Cc = 


Ma2 


ni3 


1 + c1 ■ Ma1 
c4 


nO.5 


1 + c1 ■ Ma2 






^— + c2 ■ Mai ^ 1 

h 

c3 2 

— + c2 ■ Mai 

0.5 

A2 ■ 

1 + cl ■ 


+ A1 



|Cl ) 


cl 



+ A1 ■ (1 + c1 ■ Ma1 ) 'Mai 


I 


(1 + c1 ■ Ma1 


P3 


A3 


Dc = m3 ■ v3 
0 = Ac + Be - Dc + Cc 
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APPENDIX D 


TEMPERATURE DISTORTION (EES CODE) 


m4d 

= 4.381 

Pt4d 

= 76803 

D4 = 

11 

12 ■ 

3.28 

A4 = 

K ■ 

D4^ 

4 

K = 

1.4 


R = 

287 


cp = 

p . 

K 

n 

K - 1 

Td = 

29.4 


Tt4d 

= Td 

+ 273.15 

Pta = 

: 101325 

Tta = 

Tt4d 


h = 

- (K 

: + 1 ) 

2 ■ ( 

K - 1 ) 


m4d 


Pt4d 


A4 

Tt4d ° ® 


■ K ■ 

0.5 


r /K - 1 ) 

p 

R 


Ma4d ■ 

' M 2 

■ Ma4d 




L \ J 

. 


Pt3d = Pt4d 
Tt3d = Tt4d 
18 

D3 = - 

12 • 3.28 


A3 



K 


m3d 


m3d 


m4d 


Pt3d 


A3 

Tt3d ° ® 


FK 

R 


Ma3d 



Ma3d 


h 
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a = 0.08 


A2 

= a ■ A3 

A2 

= A1 ■ 

Tt1 

= Tt4d 


1 - a 


t = Tratio 
Tt2 = t ■ Tt1 
1 - a 


L = 


Tt3 = 


t 


n = 

Pt3 

Pt4d 


■ L + 1 

Variable 

= n 


+ t 


0.5 


. 0.5 


L + 1 


Tt1 


m3 = m4d 


m3 = Pt3 


Tt4d 


nO.5 


Tt3 

A3 


Tt3 


0.5 


0.5 


m1 = t 


0.5 


■ L + 1 


Ma3 


m3 


1 + 


K - 1 'I 


m3 = m2 + ml 


Tt3 = T3 


1 + 


K - 1 


Ma3 


vS = Ma3 • ( K ■ R ■ T3 

Pt3 = P3 

K 


0.5 


K - 

1 1 

1 P 

1 + - 


• Ma3 

\ 2 

J 



K - 1 


■ Ma3 ^ 
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■ Ma1 ■ (K' R' Tt 1 )“'^ + ni2 ■ Mai ■ (K ■ R ■ Tt 2 )“'^ 


1 + 


K - 1 


0.5 


- ni3' v3 + A3 


A1 


Ttl 

R ■ — 

K 


\0.5 


■■ P3 


1 + 

K - 1 1 

■ Mal^ 


2 

, 




Ptl 


m1 

"aT 


Tt1 


nO.5 


Ma1 


1 + 


/ K - n 




nh 


Ma1 


Pti = 

P1 ■ 

r 

Ma2 = 

Ma1 


R2 = 

((K + 

1) ■ 

R2 = 

t ■ Ra 


Ra = 

((K + 

1) ■ 


K - 1 


Ma1 


K - 1 


Ma2 


1 + K ■ Ma2 ) 


Maa 


1 + K ■ Maa ) 


F2 


Fa = 


K + 

1 

r 

+ 

C\J 

Ma2 ^ '|(k - 1 )■ 

1 + K ■ 

Ma2 ^ 

. 

( K + 1 

J J 





1 ^ 1 

K + 

1 

■ 

+ 

C\J 

Maa ^ \K - 1 j] 

1 + K ■ 

Maa ^ 


1 

+ 


Pt2 = Pt1 


g = 


F2 

Fa 


q 

Tt2 = Tta + — 
cp 


Qj = q ■ m2 
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Pt2 

Pta 


j 


1 + 

■ K - 1 ■ 

• Ma2 ^ 

2 

Ma2 ■ 

(K ■ 

R 

■ T2 ) 


Tt1 


1 + 

■ K - 1 ■ 

■ Ma1 ^ 

2 

Ma1 ■ 

(K ■ 

R 

■ T1 ) 
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APPENDIX E INLET VELOCITY RATIO 


A. PRESSURE DISTORTION 


In this case, 


rtj = Tt^ 


p = p 

'*1 ^2 


Pt^ = TUPt^ 

From the definition of stagnation temperature, 


Tt —T 
F = F Ll - 2 

Tt,- T, 


Since 


( P. \ r ( P \ r 

T 2 =Tt 2 —— , and T,=Tt, — 


using equations (1), (2), (3) and (5) in equation (4), 


V,=V, - -^ 1-Ft +1 


Then, from Equation (6), the minimum value of n (when V 2 = 0) is given by 




B. TEMPERATURE DISTORTION 


In this case, 


Pt, = Pt^ 


51 



(9) 


Tt^ = fTt^ 

and 

P = P 

Since, from equations (8) and (10), 

Pti _ Pt2 
P, P2 

then Ml = M 2 and 

T T 

Re-writing equation (4) as 


( 10 ) 


( 11 ) 


( 12 ) 



and using equation (12), then 


(13) 


(14) 
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